Break up of the heavy electron at a quantum critical point 
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The point at absolute zero where matter becomes unstable to new forms of or- 
der is called a quantum critical point (QCP). The quantum fluctuations between 
order and disorder 1-5 that develop at this point induce profound transformations 
in the finite temperature electronic properties of the material. Magnetic fields 
are ideal for tuning a material as close as possible to a QCP, where the most in- 
tense effects of criticality can be studied. A previous study 6 on theheavy-electron 
material YbRh 2 Si2 found that near a field-induced quantum critical point elec- 
trons move ever more slowly and scatter off one-another with ever increasing 
probability, as indicated by a divergence to infinity of the electron effective mass 
and cross-section. These studies could not shed light on whether these proper- 
ties were an artifact of the applied field'' 8 , or a more general feature of field-free 
QCPs. Here we report that when Germanium-doped YbRh 2 Si2 is tuned away 
from a chemically induced quantum critical point by magnetic fields there is a 
universal behavior in the temperature dependence of the specific heat and re- 
sistivity: the characteristic kinetic energy of electrons is directly proportional 
to the strength of the applied field. We infer that all ballistic motion of elec- 
trons vanishes at a QCP, forming a new class of conductor in which individual 



electrons decay into collective current carrying motions of the electron fluid. 

Recent work 6 on the heavy electron material YbRh 2 Si 2 9 has demonstrated that a mag- 
netic field can be used to probe the heavy electron quantum critical point. This material 
exhibits a small antiferromagnetic (AFM) ordering temperature T N = 70 mK (Fig. la) 
that is driven to zero by a critical magnetic field B c = 0.66 T (if the field is applied 
parallel to the crystallographic c-axis, perpendicular to the easy magnetic plane) 6 . For 
B > B c , a field-induced Landau Fermi Liquid (LFL) state characterized by Ap = AT 2 
(where Ap{T) = p{T) — p Q is the temperature dependent part of the electrical resistivity) is 
established below some cross-over temperature T (B) which grows linearly with field. The A 
coefficient, being proportional to the quasiparticle-quasiparticle scattering cross section was 
found to diverge as A(B) oc 1/(B — B c ) for B — > B c . Comparative studies of the resistivity 
and the electronic specific heat C e \(T) = 7 (£?)T in the field ranges 0.5 T - 4 T (B _L c with 
B c = 0.06 T) and 2 T - 6 T (B \\ c) revealed a field-independent ratio A/^y 2 slightly smaller 
than the empirical Kadowaki- Woods ratio 13 that holds for LFL systems. This seemed to 
suggest a divergence of the effective quasiparticle (QP) mass as 1/(B - B c ) 1/2 as B -> B c . 
In this letter, we report the first-ever observation of the divergence of the QP mass at a 
QCP, established very close to B = 0. 

By alloying YbRli2Si2 with Germanium, using a nominal concentration x = 0.05, we 
have been able to fine-tune, in a new set of studies, the Neel temperature of this material 
and the critical field far closer to zero, to a point where, for the first time, we may reliably 
probe the zero-field transition using field-tuning. The phase diagram for a high-quality 
YbRh 2 (Sio.95Geo.o5)2 single crystal is shown in Fig. lb. NFL behavior dominates over a 
funnel-shaped region of the T-B phase diagram down to the lowest accessible temperature 
of 20 mK. The critical field has been suppressed to as low as B c = 0.027 T (B _L c). As in the 
undoped material, there is a broad cross-over regime between the NFL and field polarized 
LFL regime with a mean cross-over temperature T that is seen to rise linearly with the field 
B. Very weak AFM order develops in the x = 0.05 sample below Tn = 20 mK, as evidenced 
by the extremely weak anomaly in the electronic specific heat coefficient (Fig. 2a). 
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Past experience 7,8 suggested that a finite field quantum critical point has properties 
which are qualitatively different to a zero field transition, shedding doubt on the reliability 
of these measurements as an indicator of the physics of a quantum phase transition at zero 
field. However, the zero-field properties of YbRh 2 (Sii_ x Ge x ) 2 above T ss 70 mK for the 
undoped (x = 0) and doped (x = 0.05) crystals are essentially identical (Fig. 2a), suggesting 
that by suppressing the critical field we are still probing the same quantum critical point. 
In both compounds, the ac-susceptibility follows a temperature dependence x -1 ^ T a from 
0.3 K to < T < 1.5 K, with a = 0.75 14 , and the coefficient of the electronic specific heat, 
C C \(T)/T, exhibits 9 a logarithmic divergence between 0.3 K and 10 K. However, in the 
low-T paramagnetic regime, i. e. , T N < T < 0.3 K, the ac-susceptibility follows a Curie- 
Weiss law (inset of Fig. 2a) with a Weiss temperature Q w rs 0.3 K, and a surprisingly large 
effective moment fj, e s ~ 1.4/i B /Yb 3+ , indicating the emergence of coupled, unquenched spins 
at the quantum critical point. The electronic specific heat coefficient, C C \(T)/T, exhibits a 
pronounced upturn below 0.3 K (Fig. 2a). 

We now discuss the field dependence of the electronic specific heat in YbRh 2 (Si .95Ge .o5)2 
in more detail. In these measurements, magnetic fields were applied perpendicular to the 
crystallographic c-axis, within the easy magnetic plane (Fig. 2b). At fields above 0.1 T, C e \/T 
is weakly temperature independent, as expected in a LFL 15 . A weak maximum is observed 
in C C \(T)/T at a characteristic temperature T (B) which grows linearly with the field (inset 
of Fig. 3a), indicating that entropy is transferred from the low-temperature upturn to higher 
temperatures by the application of a field B > B c = 0.027 T. As the field is lowered the 
temperature window over which C e \(T, B)/T = 7 G (-B) is constant shrinks towards zero and 
the zero-temperature 7o(-B) diverges (Fig. 3a). For example, in a field of 0.05 T a constant 
value 7o(-B) m 1.54(7) Jmol~ 1 K -2 only develops below 40 mK. These results indicate the 
formation of a field-induced LFL state at a characteristic scale T <T Q (B). As the window of 
LFL behavior is reduced towards zero, an ever increasing component of the zero-field upturn 
in the specific heat coefficient is revealed in the temperature dependence. This confirms that 
the major part of the upturn in the specific heat coefficient observed in zero field is electronic 
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in character, and must be associated with the intrinsic specific heat at the QCP. 

This conclusion is also supported by the electrical resistivity data which reveal a field- 
dependent cross-over from a T-linear resistivity at high temperatures, to quadratic behavior 
Ap = A(B)T 2 at sufficiently low temperatures. Most importantly, the data show that at low 
fields and temperatures, the same scale T (b) oc b (where b = B — B c is the deviation from the 
critical field) governs the cross-over from LFL to NFL behavior in both the thermodynamics 
and the resistivity. This can be quantitatively demonstrated by noting that the finite field 
transport and specific heat data collapse into a single set of scaling relations (see Fig. 3 
inserts) , 



where <&(#) ~ (max(x, l)) -1 ^ 3 and F(x) ~ x/max(x,l). The NFL physics is described 
by the x — > oo (T >> T (b)) behavior of these equations, where dp/dT is constant and 
Cy/T oc T~ 1//3 . By contrast, the field-tuned LFL is described by the x — > limit of these 
equations. Were there any residual pockets of LFL behavior that were left unaffected by 
the QCP, we would expect a residual quadratic component in the resistivity, and the data 
would not collapse in the observed fashion. We are thus led to believe that the break up of 
the LFL involves the entire Fermi surface. 

From the second scaling relation in (1), we see that the A-coefficient of the T 2 term 
to the resistivity diverges roughly as 1/b, a result that is consistent with earlier measure- 
ments on pure YbRh 2 Si2 carried out further away from the QCP coefficient 7 (6) grows as 
& _1 / 3 (Fig. 2a). Notice that the field dependence at absolute zero temperature can be inter- 
changed with the temperature dependence at B — B c , but only in the upturn region. At 
high magnetic field deviations from the QCP, earlier measurements showed 6 that the Kad- 
owaki Woods ratio 13 K = A/rf is approximately constant. Closer to the QCP, where the 
scaling behavior is observed, K = Aj^ 2 ss b^ 1 ^ 3 is found to contain a weak field dependence 



We now turn to discuss the broader implications of our measurements. The observed 




(1) 



(Fig. 3b). 
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divergence of both the ^-coefficient and the coefficient 7 G of the T-linear specific heat cer- 
tainly rule out a 3D SDW scenario, which predicts that both quantities will remain finite at 
sufficiently low temperature in the approach to a zero-field QCP (B — > B c — > 0), but it can 
be used to obtain still more insight into the underlying scattering mechanisms between the 
quasiparticles. In a 2D SDW scenario, the scattering amplitude between two heavy electrons 
is severely momentum dependent. When used to compute the transport relaxation rate, the 
SDW scenario leads to the result A oc l//t 2 , with k the inverse correlation length 16 . The 
observed divergence in A{b) would require k 2 oc b. However, the fluctuations of the soft 
2D spin fluctuations only produce a weak logarithmic renormalization in the heavy electron 
density of states, measured by the specific heat coefficient, 7 Q oc In Thus the 2D SDW 

scenario predicts a weak divergence in the in-T linear specific heat, but a strongly field 
dependent enhancement of the Kadowaki Woods ratio in the approach to a QCP (b — > 0), 
given by 

Isdw oc ln(l/6) and K SDW oc \ . (2) 

oln (b) 

The strong violation of these predictions by our data, presented in Fig. 3a and 3b, rules 
out 2D spin fluctuations as the driving force behind the thermodynamics and the dominant 
source of scattering near the heavy electron QCP. 

Taking a more general view, scaling behavior of the transport scattering rate tells us 
that the only scale entering into the density of states and the scattering amplitude is the 
single scale T oc b of the heavy electron fluid. A truly field-independent Kadowaki Woods 
ratio would indicate that the quasiparticle scattering amplitude has the form 

A*=T* F F[{k in }^{k 0Ut }}. (3) 

The weak field dependence of the Kadowaki Woods ratio over a wide range of fields implies 
that the characteristic length scale of the most singular scattering amplitudes renormalizes 
more slowly in the approach to the QCP than expected in a SDW scenario. 

Our data also provide some insight into the thermodynamics in the vicinity of the QCP. 
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By integrating the scaling form (1) for the specific heat over temperature, the entropy 
S(T) = J* T dT'(C v /T') in the vicinity of the QCP can be described by the form 



where 77 = 1/3. The appearance of a field-dependent pre-factor in this equation forces 
the entropy to vanish at the QCP, as required by the third law of thermodynamics. The 
exponent in the pre-factor also determines the effective Fermi temperature T F (b) oc 7(6) _1 oc 
T (b) v . Thus the requirement that the entropy vanishes at the QCP (rj < 1) prevents a 
direct proportionality between the Fermi temperature of the heavy LFL and the scale T (b) 
governing the cross-over to NFL behavior. It follows that the Fermi temperature and cut-off 
temperature T Q (b) must obey a relationship of the form 



where T A is an upper cut-off that we might identify with the single ion Kondo temperature 
of the Yb 3+ ions(~ 25K). Such a power law renormalization of the characteristic energy 
scale would be expected in the presence of locally critical fluctuations that extend down 
from Ta to the infra-red cut-off provided, in this case, by the magnetic field 17 . 

In this respect, our results support the conclusions recently drawn from earlier measure- 
ments on the quantum critical material CeCu 6 _ x Au x (x = 0.1) 18 , and used in a recently 
proposed theory for quantum criticality by Si et al. 4 , suggesting that the most critical scat- 
tering is neither three, two or even one dimensional, but local- as if the most critical fluc- 
tuations in the underlying quantum phase transition are fundamentally "zero dimensional" 
in character. 

One of our most striking observations is that below T as 0.3 K where x{T) follows a Curie- 
Weiss law, the electronic specific heat coefficient C e \(T)/T for both samples starts to deviate 
towards larger values, separating away from the — log T dependence that is valid 9 up to 10 K 
(Fig. 2a). This "upturn" continues in the x = 0.05 sample down to approximately 20 mK, if 
the critical field of 0.027 T (B ± c) is applied. We ascribe this intrinsically electronic feature 




(4) 




(5) 
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to the critical fluctuations associated with the zero-field quantum phase transition that exists 
at a slightly larger Ge concentration. The unique temperature dependence of C e \(T)/T for 
T < 0.3 K is disparate from the linear temperature dependence of the electrical resistivity 
which holds all the way from > 10 K to T m 10 mK. Since the former (thermodynamic) 
quantity probes the dominating local 4/ ( "spin" ) part of the composite quasiparticles , while 
the latter (transport) quantity is sensitive to the itinerant conduction-electron ("charge") 
part, one may view the observed disparity as a direct manifestation of the break up of the 
composite fermion in the approach to the QCP. 
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FIG. 1. Evolution of e, the exponent in Ap(T) = [p(T) — p a \ oc T £ , within the temperature-field 
phase diagram of YbRh2(Sii - x Ge x )2 single crystals. The non-Fermi-liquid (NFL) behavior, e = 1 
(yellow color), is found to occur at the lowest temperatures right at the quantum critical point 
(QCP), B = B c [ a, x = 0, B c = 0.66 T (|| c), residual resistivity p Q = l^$7cm; b, x = 0.05, 
B c = 0.027 T (_L c), p = 5pftcm ] and in a largely extended field range at higher temperatures. 
For B > B c , a broad cross-over regime from the NFL state to the field-induced heavy Landau 
Fermi-liquid (LFL) state (at lower temperature) is stated. The LFL state is characterized by 
Ap(T) oc T e , e = 2 (blue color). As shown in (a) the antiferromagnetically ordered phase of pure 
YbRli2Si2 below Tn = 70 mK and B c shows, owing to an extremely small ordered moment, the 
outward appearance of a heavy LFL state, too. Its phase boundary to the paramagnetic state is 
manifested by a rapid change in e from 2 to 1. 

The low ordering temperature of pure YbRli2Si2 increases as external pressure is applied 9 . The 
extrapolation of T^(p) — > yields a critical pressure p c = —0.3(1) GPa, reflecting that a small 
expansion of the unit cell volume, V, would tune Tn — > 0. This can be achieved by the substitution 
of Si by the isoelectronic, but larger, Ge 10 . Studies of the electrical resistivity under pressure 
revealed a Tn oc (p + p c ) n variation, with n = 1.33 for both compounds. The T^{p) -dependences 
of the x = and x = 0.05 crystals can be matched if all x = 0.05 data points are shifted by 
the same amount Ap = —0.17(2) GPa to lower pressure 10 , yielding Tn = 20(5) mK. Using the 
bulk modulus B = 189 GPa of YbRh 2 Si 2 n , the small pressure shift of Ap = -0.17(2) GPa is 
equivalent to a volume expansion of AV = 0.14(3) A 3 . This transforms into an effective Ge content 
x cS = 0.019(6), if the value AV/V = 7.65(78) % for the relative change of the unit-cell volume with 
Ge concentration in the solid-solution YbRl^Sii-zGe^ is used, with V(x = 0) = 158.4(2) A 3 
and V(x = 1) = 166.07(54) A 3 , cf. Ref. 12 , in agreement with microprobe analysis 14 . 
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FIG. 2. Low-temperature electronic specific heat of YbRI^Sii-^Gea;^ single crystals as C c \/T 
vs T in semi-log plots at zero field and at low values of the applied magnetic field B. Insets show 
low-T B = ac-susceptibility as x~ l vs T (a) and magnetization as M vs B (b). C e \ is obtained by 
subtracting the nuclear quadrupolar contribution, Cq = ckq/T 2 (with oq = 5.68 x 1CT 6 JKmol -1 , 
calculated from recent Mossbauer results 11 ) (a) and, in addition, the nuclear Zeeman contribu- 
tion C hf = a{B)/T 2 (b), from the raw data. Here, a(B) has been deduced by plotting CT 
vs T 3 . The magnetization, M vs magnetic field B (black points in the inset), is calculated via 
(Bhf — B)/A, with A the hyperfine coupling constant for Yb in this compound and the hyperfine 



field Bhf = (a(B) — oiq) / cx.dip\ adip represents the strength of the nuclear magnetic dipolar inter- 
action and amounts to 7.58 x 10~ 8 JKmol _1 T -2 , Ref. 19 . With the assumption of A = 120 T/// B , 
the data points agree perfectly with the measured magnetization curve at 40 mK (red line in the 
inset of (b). 

The B = results shown in (a) reveal an upturn in C e \(T)/T for paramagnetic YbRl^Sii-zGe^ 
(x = 0, Tn = 70 mK; x = 0.05, Tn = 20 mK) below T = 0.3 K. In the same temperature range 
the susceptibility x{T) shows a Curie- Weiss law, \ l oc (T — 0) [inset of (a)]. For both samples 
very similar values are found for the Weiss temperature, « —0.3 K, as well as for the large 
effective moment, iJ e g 1.4^B/Yb 3+ . For YbRh2(Sio.95Geo.o5)2, entropy is shifted from low to 
higher temperatures when a magnetic field is applied (b). The cross-over temperature between 
the field-induced LFL state (C e \(T)/T ss const.) and the NFL state at higher temperature is de- 
picted by the position of the broad hump in C e \(T)/T which shifts upwards linearly with the field, 
B < 0.8 T (_L c). 
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FIG. 3. Field dependences of the Sommerfeld coefficient, j , of the electronic specific heat 
(a) and of the ratio of the ^-coefficient in the T 2 term of the electrical resistivity and 7 2 (b) 
for YbRh2(Sio.95Sio.o5)2- Note that 7 D and A are proportional to the effective quasiparticle mass 
and the effective quasiparticle-quasiparticle scattering cross section, respectively. The magnetic 
field was applied perpendicular to the c-axis, and the applied field values are corrected, on the 
abscissae, by the value of the critical field, B c = 0.027 T. 7 -values in (a) were obtained from two 
different samples: Three independent measurements on sample #1 are displayed by closed symbols 
(circles, up and down triangles). The open symbols (diamonds) show the results of measurements 
on sample #2. As B — > B Cl 7 D diverges oc (B — 5 C )~°' 33 (red line), i.e., much stronger than 
logarithmically (black dashed line). The symbols used in the semi -log plot K = Aj-fl vs (B - B c ) 
of (b) correspond to values for the electronic specific heat coefficient shown in (a). Half filled 
circles (squares) display data for which the ^-coefficient of the electrical resistivity was determined 
by extrapolating (interpolating) A(B — B c ) with respect to (B — B c ). The half filled diamond 
represents a point for which the 7 G value was obtained by interpolation. The black dashed line 
indicates K S dw oc {B - B c ) ln 2 (B - B c ) ~ for the 2D SDW scenario 6 . This is at strong variance 
from the (at B — B c < 0.3 T) experimentally observed K oc (B — B c )^ 1 ^ , arising from the 
stronger than logarithmic increase of 7 D upon cooling (red line in a). For (B — B c ) > 0.3 T, 
K becomes field independent within the error bars at a constant value of 5.4 ^f]cmmol 2 K 2 J~ 2 
(green horizontal line in b). A similar high field behavior has been reported previously 6 on pure 
YbRli2Si2. Noteworthy, an almost identical value for K was found. Insets show scaling behavior 
of the low-T electronic specific heat where, according to equation (1), the ordinate is displaying 
Q(B, T) = (B — i? c ) ' 33 C e \/T, a, as well as of the temperature derivative of the electrical resistivity, 
dp/dT, b, as a function of T/(B - B c ). 
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